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Abstract. We use the semi-analytic method developed by Fedeli et al. for computing strong-lensing optical depths to study 
the statistics of gravitational arcs in four dark-energy cosmologies. Specifically, we focus on models with early dark energy 
and compare them to more conventional models. Merger trees are constructed for the cluster population because strong cluster 
lensing is amplified by factors of two to three during mergers. We find that the optical depth for gravitational arcs in the early 
dark-energy models is increased by up to a factor of ~ 3 compared to the other models because of the modified dynamics of 
cluster formation. In particular, the probability for gravitational arcs in high-redshift clusters is considerably increased, which 
may offer an explanation for the unexpectedly high lensing efficiency of distant clusters. 



1. Introduction 

While the presen t domi n ance of dark ene r gy is well es tablished 
JGoldstein et alJ I20031 lHawkins et all I20031 ISnergel et alJ 



2003I I2004 iReholoetalJ I2004 iReadhead et alJ |2£ 
Riess etal.ll2004l iTe gmarket al.ll2004h . its evolution in par- 
ticular in the early universe is largely unconstrained. An 
interesting class of models for dynamical dark energy is 



characterised by 

energy a t early times ( i Ferreira & Jo^ce 
l2001alht ICaldwell et alJ l2003t IWettericr 



low but non- vanishing density of dark 
" 19981 iDoran et alJ 



2004). Non-linear 



structure for mation has recently been studied for this class 
of models by iBartelmann et alJ J2005t) in the framework of 
the spherical collapse model. Interestingly, it was found that 
non-linear structures are expected to form substantially earlier 
in such early dark-energy models if they are normalised so as 
to be compatible with the large-scale temperature fluctuation 
amplitude of the cosmic microwave background. For two 
specific models with early dark energy, the population of 
galaxy clusters is expected to evolve by approximately an 
order of magnitude less strongly than in the standard ACDM 
model. 

Should this come close to reality, a rich population of mas- 
sive galaxy clusters would be present at high redshift which 
is completely unexpected in ACDM. Similarly, the dynamical 
activity within the cluster population due to substantial merg- 
ers with sub-halos would be shifted or extended towards higher 
redshift. 

The problem of the non-linear evolution of cosmic struc- 
tures in presence of dark energy has been recently addressed 
also from a more general point of view b y several authors. Fo r 
instance. iMota & van de Bruckl |2004l) IZeng & Gaol j2005), 
iMaor & Lahavl J2005t) and Wan j J2006I) analyse different as- 
pects of this issue for both constant and time dependent dark 
energy equation of state parameter, allowing for dark energy 



clustering and coupling to dark matter. They outline very dif- 
ferent properties of the virialized objects dependi ng on the be- 
haviou r of the dark energy fluid . Additionally, IZeng & Gad 
J2005t) and lManera & Motal J2005I) explore the outcome of this 
different non-linear evolution on the predicted number counts 
for high mass dark matter haloes (galaxy clusters), finding 
several significant effects. They discover in particular that the 
number counts of massive structures increase if small scale 
clustering of dark energy is allowed, while it decrease if the 
amount of dark matter coupled to dark energy grows. 

One interesting and due to its non-linearity highly sensi- 
tive way for probing the massive end of the cluster popula- 
tion is the st rong lensing effect. Although the issue is still 
contro v ersial (IBartelmann et al]ll998t IMeneghetti etaDfe oOO , 
1 20034 IWamhsganss etalJl2004 lOalal etalJ 1200.4 iFi er.alJ 
2005t iHennawi et alJ 12005 i. it seems to be at least difficult 
within the ACDM model to reproduce the observed abundance 
of strong-lensing events in cluster cores, the s o-called gravi- 
tation a l arcs. Arcs in clusters a t high redshift jHas inger et all 
199 81 iThompson et alJ 1200 it IZaritskv & Gonzale j 120031 
iGladders etall2 003) are similarly puzzling because they indi- 
cate that even clusters at z > 1 can already be concentrated and 
massive enough to be strong gravitational lenses for a source 
population that is not too distant from them. 

Dynamical activity in galaxy clusters was identi- 
fied as highly important for their strong-lensing abil ities 
( Bar telmann et alJll99^: IMeneghetti eTat] l2003bt iTorri et alJ 
20041- lFedeli et alJ2 006). The enhancement of the gravitational 
tidal (shear) field while clusters are merging with massive 
halos can transiently, but substantially increase their strong- 
lensing cross sections. As much as about half of the total 
optical depth for strong cluster lensing may be contributed by 
merging clusters. The effect is strong because mergers can turn 
clusters into strong lenses that would otherwise fall below the 
threshold because they are not massive or compact enough. 
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Major cluster mergers thus open the huge, exponentially rising 
reservoir of moderately massive clusters for strong lensing. 

Cosmological models reconciling an appreciable cluster 
abundance at high redshift, and thus also a high level of dy- 
namical cluster activity, with independent cosmological con- 
straints e.g. from the CMB are thus particularly interesting in 
view of strong cluster lensing. Sufficiently detailed numerical 
simulatio ns are costly and beyond scope for a parameter study. 
Recently, Fed eli et ail J2OO6) have developed a (semi-)analytic 
method for computing strong-lensing cross sections for galaxy 
clusters with and without taking cluster mergers into account. 
This method opens the way to systematically test a variety of 
cosmological models for their consequences for strong cluster 
lensing. We use it in this paper to study the statistics of strong 
cluster lensing in two exemplary cosmologies with early dark 
energy and compare them to the ACDM model and a model 
with a constant equation-of-state parameter w > - 1 . We do not 
focus on other aspects related to the dark energy fluid behaviour 
like small scale clustering or coupling to dark matter. 

Section 2 reviews extended Press-Schechter theory as it 
will be needed later, and Sect. 3 summarises the cosmological 
models used. The construction of merger trees and the compu- 
tation of strong-lensing cross sections are described in Sects. 4 
and 5. Section 6 outlines expectations, Sect. 7 quantifies the 
results, and Sect. 8 summarises and concludes the paper. 

2. Press-Schechter theory 

We begin by reviewing the basic features of the excursion- 
set approach to t he de rivation of the Press-Schechter 
JPress & Schechterl 1 1974b mass function and the con- 
ditional mass function of virialised dark- matter halos . 
Comprehensive treatments can be found in iBond et al J ill 99 lb 
and lLacev & Cold dl993h . The central physical quantity is the 
primordial Gaussian density-fluctuation field S(x), filtered on 
a scale R corresponding to a mass M. The filtered field 6m{x) 
remains a Gaussian random field whose variance S — S (M) is 
a monotonically decreasing function of mass. If the filter is a 
top-hat function in Fourier space, Sm(x) at a fixed location x 
performs a random walk as a function of M, and thus of scale 
R. 

When 6m(x) rises above a critical, redshift-dependent 
threshold <5 c (z), a halo of mass M is expected to form at the 
location x at redshift z. The threshold S c (z) is usually obtained 
from the spherical collapse model by linearly extrapolating the 
initial overdensity to the time when the collapsing halo reaches 
virial equilibrium. The probl em of a random walk w ith a fixed 
absorbing barrier is solved in lCha ndrasekharl d 1 943l) and leads, 
under the previous assumptions, to the mass function 



n(M, z) = — - 



Sc(z) 



M V2^5 3 / 2 D + (z) 



AS 



AM 



exp 



2SD + {z) 2 



(1) 



The mass function is defined such as n(M, z)dM is the comov- 
ing number density of structures with mass between M and 
M + AM at redshift z, where p m ,o is the mean matter density 
at present time. Unlike common practice, we explicitely intro- 
duced the linear growth factor D+(z) instead of incorporating 



Table 1. Cosmological parameters of the four models used in 
the present work. 





EDE1 


EDE2 


Wde = -0.8 


ACDM 




0.82 


0.78 


0.80 


0.84 


h 


0.67 


0.62 


0.65 


0.65 


£2m,0 


0.33 


0.36 


0.30 


0.30 




0.67 


0.64 


0.70 


0.70 



it into the critical overdensity in order to emphasise that the 
redshift dependence in 5 c {z) is exclusively due to the evolution 
of the spherical collapse model with redshift. In an Einstein-de 
Sitter universe, 5 C ~ 1.686 is a constant, and it evolves only 
gently in a ACDM universe, but it changes considerably in par- 
ticular in the cosmological models with early dark energy. We 
thus have to account for its redshift evolution. 

The formalism sketched so far can easily be extended to 
construct the conditional mass function which quantifies the 
probability for a halo of a given mass Mq at a given redshift 
z to have a progenitor of a lower mass M p at a higher redshift 
z + Az. Since the variance of the density field 5m{z) filtered on 
a scale corresponding to a given mass M decreases monotoni- 
cally with M, this is equivalent to the probability that a halo of 
variance S (Mo) at a given redshift had a higher varia nce S (M v ) 
at a h igher redshift. This probability is given by jLacev & Cold 
Il993h 



K(AS,A(o) = 



1 



All) 



V2^ AS 3 / 2 



exp 



Aco 1 



2AS 



(2) 



where AS = S(M p ) - S (Mo), and Aco represents the redshift 
step considered, 



Aco 



S c (z + Az) S c (z) 



D + (z + Az) D + (z) 



(3) 



In other words, Eq. $2^ gives the probability for a dark-matter 
halo to undergo a change in variance AS due to hierarchical 
accretion in the redshift interval Az. 

If we want the probability for the halo of mass Mq to have 
a progenitor corresponding to a change in variance lower than 
AS within the same redshift interval, we simply have to inte- 
grate the above equation, obtaining the cumulative probability 
distribution 



f / Aco \ 

/(AS, Aw) = K(Ai", Ato)AA£ = erfc -— = , 
Jo IV2AS7 



where 



erfc(x) 



2 p 



At 



(4) 



(5) 



is the complementary error function. Equation (0} is just the 
probability for the mass of the progenitor M p to be larger than 
the mass corresponding to the variance AS - S {Mq). 
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3. Dark-energy models 

Dark energy generalises Einstein's cosmological constant, re- 
placing it by a term varying with redshift in Friedmann's 
equation. Viable model universes with dark energy must be 
adapted to comply with present-day observational data. Many 
models motivated by elementary-particle physics introduce 



see for instance IWetterichll 19881 


Peebles & Ratra 1988, 2002 


IBrax & Martini 


2000), whose pressure and energy density are 



related by the perfect-fluid equation of state 



de = W de p de C 



(6) 



The parameter Wd e is typically a function of cosmic time or 
redshift. A cosmological constant has w de = — 1 at all red- 
shifts. Different models give rise to different (often not ana- 
lytic) shapes for the function w^iz), and there are various ways 
to parameterise a chosen quintessence model. 

Negative pressure at all times implies that the dark-energy 
density parameter Qd e (z) will fall to zero for increasing red- 
shift. If, however, the equation-of-state parameter w is al- 
lowed to rise above zero, models can be constructed in which 
^de(z) settles at a small positive value. The presence of a non- 
vanishing dark-energy contribution in early epochs of the cos- 
mic evolution can have many interesting consequences on the 
CMB temperature fluctuations, the geometry and the age of the 
Universe, and the linear and non-linear aspects of structure for- 
mation. Dark-matt er halos on all mass s cales ma y form sub- 
stanti ally earlier dBartelmann et alj|2005l see also iDolag et alJ 
2004), potentially causi ng large effects on the statistics of 
str ong-lensing event s (cf. Ba rtelmann et all2 003). 

IWettericff J2004) proposed a useful characterisation of cos- 
mological models with early dark energy which makes use of 
only three parameters; the present density parameter in the dark 
energy Qde,o> the present equation-of-state parameter w^ q, and 
an average value for the dark energy parameter at early (struc- 
ture formation) times, 



de.sf 



= -(lnOea)" 



J In cu 



^de(^) din a 



(7) 



where a eq is the scale factor at matter-radiation equality. 

For sufficiently low Wd e ,o, such phenomenological early 
dark-energy models reproduce the accelerated cosmic expan- 
sion in the present-day Universe similar to cosmological- 
constant models and can thus be arranged to agree with low- 
redshift observations. If Qd e , s f is small enough, they can also 
reproduce the CMB temperature fluctuations. We shall in- 
vestigate here the same two early dark-energy models as in 
iBartelmann et alJ J2005I) . which have fi de , s f = 0.04, spectral 
indices for the primordial density-fluctuation power spectrum 
of n — 1.05 (model EDE1 henceforth) and n = 0.99 (hereafter 
model EDE2). For comparison, we also include a model with 
a constant equation-of-state parameter w = -0.8 and a conven- 
tional ACDM model for reference. 

The values of the other cosmological parameters, such as 
the present-day matter-density parameter Q. m fi, the dimension- 
less Hubble constant h and the normalisation of the power 
spectrum expressed by cr 8 , were determined such as to match 



N 

QJ 

■a 



-0.4 



-0.6 



-0.8 




_l I I I I I I I I I I L 



_L 



_1 I I I I I I L 







2 3 
z 



Fig. 1. The redshift evolution of the equation-of-state parame- 
ter for the four cosmological models used in this paper. These 
are a ACDM model (green dashed-dotted line), a model with 
constant w^ e = -0.8 (blue long-dashed line), and two early 
dark-energy models with different spectral indices for the pri- 
mordial density fluctuations (the black solid line represents 
EDE1 and the red-dashed line EDE2). 



the power spectrum of th e CMB temperature fluctuations 
(Soer gel et all I2003L 120061), constraints from th e large-scale 
structure of the Universe jTegmark et alJ EoO^. and obser- 
vations of type-la supernovae (IRiess e t al. 2004). The values 
characterising the four cosmological models used in this paper 
are listed in Tab.^ while Fig^shows the redshift evolution of 
the equation-of-state parameter in these cosmologies. 

As can be noted, the early dark-energy models approach 
close the cosmological-constant scenario at very low redshifts. 



4. Merger trees 

We now proceed to use the extended Press-Schechter formal- 
ism summarised in Sect. 2 for a Monte-Carlo realisation of 
merger trees. The procedure is qu ite straightforward, and we 
refer to Somerville & Kolatt (1999) for a detailed discussion 
and to Randal l et alJj2002l) and lCassano & Brunettil J2005I) for 
some applications. 

4.1. Monte-Carlo simulations 

Consider a halo of mass Mq at the present time (z = 0). If we 
draw a random number r in the interval [0, 1] and solve the 
equation J(AS,Acj) - r — 0, we draw a value for the change 
in the variance corresponding to the halo compliant with the 
merger rate ©. 
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Given the variance 5(Mq) of the halo's original mass, we 
obtain a new value of the variance and convert it to a new 
mass which is the mass of the progenitor M p . If we choose 
a sufficiently small time interval, we can assume that the entire 
change of the halo's mass is due to a unique, binary merging 
process with another halo of mass AM = Mo - M p . If we re- 
peat this process for earlier progenitors at subsequent redshift 
steps, we obtain the merger history of the original halo up to a 
given redshift. At the end of this procedure, we have obtained 
the value of the halo's mass and that of its progenitors for each 
redshift step, i.e. a merger tree. 

The choice of the time interval needs some care. It has to 
be small to justify the assumption of binary mergers, but not 
too small to avoid that the results be dominat ed by numeri- 
cal no ise. Following the rule-of-thumb given byE acev & Cold 
il 19931) . we use a time step such that 

(see also ISomerville & Kolattl fl 999t) where AM C is the mass 
of the smallest sub-halo required to be resolved individually. 
If M p or AM fall below AM C , the process does not represent 
an individual merger, but smooth accretion. It follows from the 
above expression that the lower initial masses Mo require larger 
time steps. 

A set of Monte-Carlo realisations of merger trees is suc- 
cessful if the population of structures that it produces agrees 
with the theoretical mass function at any given redshift. As 
ISomerville & Kolattl {T999) pointed out, this is not strictly so 
if we consider only binary mergers and smooth ac cretion as 
we are doing here. Several authors ( Ben son et alJ 120051) ar- 
gued that this may be due to an intrinsic inconsistency in the 
exten ded Press-Schechter formalism, and ISomerville & K olatt 
(1999) suggest that the problem can be mitigated considering 
multiple mergers and smooth accretion. Nonetheless, the differ- 
ence between the halo-mass distributions following from the 
merger-tree simulations and expected from the mass function 
is significant only at redshifts beyond our interest, and we con- 
firmed with several tests the good agreement between the two 
halo-mass distributions. 

4.2. Our sample 

We consider a sample of N — 500 dark-matter halos whose 
present-day masses are uniformly distributed within M m f 
!() '■// XI and M sup = 2.5 x 1O 15 A _1 M . It is plausible that 
structures with mass below Mm at z = do not contribute ap- 
preciably to the to tal lensing efficiency (see the discussion in 
iFedeli et all2006l) . For each halo, we compute the appropriate 
time step from (|8j and split it into two progenitor halos. Then, 
we proceed with the more massive progenitor as the starting 
point for the next step. We repeat this procedure until the red- 
shift exceeds the source redshift z s (which is chosen individu- 
ally for each halo in the sample, see Sect. 4 for details) or the 
mass of the halo falls below AM C . 

We show in Fig.|2]the merger histories (that is the evolution 
of mass with redshift) of five halos selected from our sample of 
500 halos for a ACDM model. Sudden discontinuities in the 




0.5 1 1.5 2 



z 

Fig. 2. Exemplary merger histories for five dark-matter halos 
randomly selected from our sample in a ACDM universe. The 
merger histories are extended up to the source redshift for each 
individual halo. 

mass are evident, each of which corresponds to a merger be- 
tween the main halo and a massive sub-halo. 

5. Strong-lensing statistics 

In order to compute the efficiency of dark-matter halos as 
strong cluster lenses, specifically for producing long and thin 
arcs, we model each halo as a NFW density profile with ellip- 
tically distorted lensing potential. Following Meneg hetti et alJ 
( 20033), we adopt an ellipticity for the iso-potential contours 
equal to e = 0.3 for all halos. Deflection-ang le maps for 
such a lens model can be calculated analytically (Bartelmann 
199 8). We then use th e fast, semi-analytic method developed 
bv lFedeli et ail feOOfj) to compute the cross sections. We de- 
scribe here only its main features and refer the reader to the 
cited paper for details. 
The lens equation 

y(x) =x- a(x) (9) 

relates the (dimensionless) position y of a point source to the 
positions x of its images. In the single lens-plane case we con- 
sider here the Jacobian matrix of the mapping (J9jl is symmetric, 
hence it can be diagonalised through an orthogonal transforma- 
tion. The eigenvalues of the Jacobian give the distortions of the 
image along the independent directions of the corresponding 
eigenvectors. Thus, the ratio of the eigenvalues determines the 
length-to-width ratio of the images of point-like sources. 

Based upon this consideration, the cross section for gravita- 
tional arcs with length-to-width ratio exceeding a given thresh- 
old can be evaluated as the integral of the inverse of the mag- 
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nification (the Jacobian determinant) over the region of the 
lens plane where the ratio of the eigenvalues is larger than 
the threshold. The region of integration will of course be a 
stripe surrounding the critical curves of the mapping. However, 
this calculation ignores the fact that real sources are extended 
and non-circularly shaped, which both enlarge the cross section 
compared to point sources. 

Thus, to obtain the length-to-width ratio of extended 
sources, we convolve the ratio of the eigenvalues with a step 
function of width equal to the source size (assumed to be 
0. 5" in radius), and we use the elegant formalism developed 
in K eetonl J200ll) to account for source ellipticities. Source el- 
lipticities are randomly drawn from a flat distribution between 
0.5 and 1. 

This method allows us to calculate the cross section <x d for 
arcs with length-to-width ratio exceeding some given threshold 
d. We choose d — 7.5 here and show one plot with d — 10 for 
comparison later. 

These semi-analytic cross sections are in excellent agree- 
ment with the results from fully numerical ray-tracing simula- 
tions. Moreover, their computation is substantially faster since 
the method does not require costly operations such as finding 
all images of every source and refining the source distribution 
near caustics on an adaptive grid. 

We calculate cross sections both ignoring and accounting 
for merger processes which transiently increase the lensing ef- 
ficiency. When a merger with a sub-halo of mass larger than 
5% of the main halo's mass occurs, we model the interaction 
as follows. The two clumps of dark matter approach each other 
at a constant speed starting from an initial distance set to the 
sum of their virial radii, r v 1 + r v 2. The process concludes when 
the profiles overlap completely, i.e. when their centres coincide 
in projection, and its duration is set to the dynamical timescale 



yn 



0v,i + r v ,2) 3 
\ G(Mi + Mi) 



(10) 



For a fixed and constant source redshift z s , we can compute the 
optical depth Td(z s ) once the strong-lensing cross sections for 
each halo at all redshift steps between the observer (z = 0) and 
the sources (z = z s ) are known. It is 



Td(Zs) 



= — r r 

4nD 2 s Jo Jo 



N(M, z) o- d (M, z, z s )dMdz , (11) 



where D s is the angular-diameter distance to the source sphere 
and N(M,z)dz is the number of structures with mass M in- 
cluded in the shell between redshift z and z + dz. However, in 
order to account for the source redshift distribution, we ran- 
domly assigned to each dark-matter halo an individual source 
redshift z v , i = 1,...,W, draw n from the r edshift distribu- 
tion of faint blue galaxies given in lSmail et alJ (fl995) (see also 
iBartelmann & SchneideJ200lh 



P(Zs) 



P 



zl exp 



(12) 



The parameters zo and /3 define the average redshift and the 
steepness of the high-redshift tail of the distribution, respec- 
tively. In this work, we used the conventional values zo = 1 and 



0.8 



0.6 



0.4 



0.2 
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1 1 1 1 I 




Fig. 3. The differential (black solid line) and cumulative 
(red dash-dotted line) source-redshift distributions given by 
Eq. G3- 



j8 = 3/2. Given this choice, the distribution peaks at z m — 1.21. 
Figure [3] shows this distribution together with its cumulative 
function defined by 



P(Zs)= P 
Jo 



p(z)dz . 



(13) 



Using the distribution (I12> . we can define the average optical 
depth as 



Tdfo) p(zs)dz s 





Jr»oo r /"*Zs /"*' 
Jo Jo 



cr d (M,z,z,)N(M,z) 



dMdz 



AnDi 



p(z s )dz s .(14) 



Since each halo in our study is characterised by a source red- 
shift randomly drawn from the distribution d!2i . we can omit 
the weighting with p(z$) when we discretise the integral over 
source redshift in (I14t . However, this is not possible for the 
mass integration, since the masses of the halos are randomly 
drawn from a uniform distribution, which requires the weight- 
ing with the halo mass function is necessary. 

The source-redshift distribution p{z s ) formally extends to 
an infinite source redshift, but obviously this is not true in re- 
ality. We set the maximum source redshift to z max = 7.5. As 
Fig- ID shows, the probability to find a source at this redshift 
can safely be neglected. Since we operate on a discrete sample 
of N halos, each of which is characterised by a mass M, and a 
source redshift z s ( , we can rewrite J14i as 



T d= )j JV(M,z)dM 

Jo [jrf 4nD 2 J Ml 



dz . (15) 
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13 13.5 14 14.5 15 15.5 
Log(M[M h" 1 ]) 

Fig. 4. The mass function for dark-matter halos in the mass 
range [10°, 2.5 x 10 15 ] hT l M Q at redshift z = 0.5 for the four 
cosmological models used in this paper, as labelled in the plot. 

The integrand of this equation is the optical depth per unit red- 
shift, i.e. the contribution to the optical depth from halos at dif- 
ferent redshifts, accounting for the source-redshift distribution 

V CTd(Mi,Z,Z s ,i) f ..... .... 

t d(z)=) l a — 7^2 — I N(M,z)dM . (16) 

This will be the central quantity in our strong-lensing analysis. 

6. Expectations 

Before turning to the results, it is useful to evaluate the expec- 
tations in or der to gain a better unde rstanding of the problem. 
As shown by Bartelman n et alJ (|2005), the formation of nonlin- 
ear cosmic structures occupies a larger redshift range in early 
dark-energy cosmological models. Structures form earlier and 
the formation process lasts longer. This increases the merger 
probability for a given halo at high redshift as well as the to- 
tal number of structures of a given mass which are found at a 
given redshift. Figure|4]shows the mass function Q at a fixed 
redshift z = 0.5 for the four cosmological models used in this 
paper. 

Evidently, the mass function is lowest for a ACDM model, 
and only slightly higher for a model with constant equation-of- 
state parameter w = -0.8. It is highest (by up to an order of 
magnitude at the high-mass tail) for the two early dark-energy 
models. This reflects the different halo-formation histories in 
different cosmologies. In the EDE1 and EDE2 models, struc- 
ture formation begins earlier, hence at a given (suitably low) 
redshift, the abundance of halos is larger. 

In Fig. |5] we show the merger rate for a halo of mass 
M = \Q lA h- x M Q and a sub-halo of mass M/2 as a function 



I i i i i I i i i i I i i i i | i i i i | i i i i | 




_l i i i i I i i i i I i i i i I i i i i I i i i i I 

1 2 3 4 5 
z 

Fig. 5. The probability for a dark-matter halo of mass M = 
10 14 /t'Mq to merge with a sub-halo of mass M/2 is shown 
as a function of redshift per unit logarithmic mass of the merg- 
ing sub-halo and per unit logarithmic cosmic time. Results are 
shown for all four cosmological models considered here, as la- 
belled in the plot. 

of redshift. By merger rate, we mean the probability for a halo 
of mass M to merge with a sub-halo of mass M/2 at redshift 
z per unit logarithmic sub-halo mass and per unit logarithmic 
cosmic time. It can be obtained as the appropriate limit of the 
conditional probability distribution Eq. 0- 

Regarding the merger rate, we also note the difference be- 
tween the behaviour of early dark-energy models and of models 
with constant equation-of-state parameter. At high redshift, the 
early dark-energy merger rate is significantly higher than for 
the other two models, but becomes essentially the same below 
redshift ~ 1.2. 

This can again be understood in terms of the different dy- 
namics of structure formation. Keeping the mass of the sub- 
halo fixed, we expect more halos of such mass to be available 
at high redshift with which the main halo can merge, because 
structure formation begins earlier in early dark-energy mod- 
els. On the other hand, structure growth begins later in models 
with a constant equation-of-state parameter and proceeds more 
rapidly. Thus, at a sufficiently low redshift, the abundance of 
such halos equals that in early dark-energy models, giving rise 
to an almost identical merger rate. 

It is worth emphasising here that the differences shown be- 
tween the different cosmological models are also due, in part 
or mainly, to the different normalisation <x 8 of the power spec- 
trum, which is chosen to make the models agree with the CMB 
observations. 

Recalling that the source-redshift distribution peaks at red- 
shift ~ 1 .2, we expect the different merger rates to have little 
influence on the optical depth. On the other hand, since the 
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Fig. 6. The logarithm of the optical depth per unit redshift for arcs with length-to-width ratio exceeding d — 7.5 for each of the 
cosmological models studied here (top and bottom-left panels), and for d — 10 in the model EDE1 (bottom-right panel). Black 
curves show optical depths obtained including halo mergers with sub-halos, while red curves are obtained ignoring the effect of 
halo interactions. 



optical depth is essentially an average of the cross section of 
different halos weighted by their relative abundances, we ex- 
pect the difference in the mass function to severely affect the 
strong-lensing statistics. In early dark-energy models, the opti- 
cal depth per unit redshift should exceed those in the ACDM 
model and the model with a constant equation-of-state param- 
eter of w — -0.8. 

We show in the next section how well this expectation is 
satisfied. 

7. Results 

We discuss now the expected behaviour of the optical depth 
per unit redshift (see definition in Sect. 4) in the different dark 
energy cosmologies considered in this work. The occurrence 
of gravitational arcs is highly sensitive to the abundance and 
internal structure of galaxy clusters, which in turn depends on 
the linear and non-linear evolution of density fluctuations. We 
thus expect that the presence and behaviour of dark energy can 
affect it. 

We show in Fig.[6]the optical depth per unit redshift for arcs 
with a length-to- width ratio exceeding d = 7.5, obtained for 
each of our four cosmological models. For the model EDE1, we 
also show the result for d = 10. The optical depths accounting 
for and ignoring halo mergers are compared. As expected, the 
lensing efficiency vanishes near the observer and approaching 



the source redshift because of the geometrical drop in lensing 
efficiency. 

Cluster mergers increase the optical depth per unit redshift, 
and thus also the total optical depth, factors up to 2 or 3 in all 
dark-energy models. The enha ncement due to mer gers appears 
more uniform than obtained bv lFedeli et al . (2006). This is due 
to the more than one order of magnitude larger sample used 
here and to the much higher time resolution adopted (up to 10~ 2 
in redshift). 

Quite obviously, increasing the length-to-width threshold 
decreases the lensing efficiency, but the features due to merger 
processes remain qualitatively the same. 

The main result is that mergers enhance the lensing effi- 
ciency by about the same amount for each model because the 
merger rate is almost the same in the redshift range relevant for 
strong cluster lensing. However, note that the absolute value 
of the optical depth per unit redshift is higher in early dark- 
energy models, which is better seen in Figs. and [S] There, 
we compare the optical depth per unit redshift for arcs with 
length-to-width ratios exceeding d = 7.5 in the four cosmolo- 
gies, accounting for and ignoring cluster mergers, respectively. 

This effect was also expected because of the difference in 
the abundance of halos of a given mass in various cosmological 
models. These figures show that, both with and without the ef- 
fect of halo mergers, the lensing optical depth per unit redshift 
is higher by factors up to ~ 3 in early dark-energy models com- 
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Fig. 7. The logarithm of the optical depth per unit redshift for 
arcs with length-to-width ratio exceeding d — 7.5 obtained for 
the four dark-energy models considered here. The lensing ef- 
ficiency shown here takes transient boosts by cluster mergers 
into account. 

pared to the other models. At redshifts above ~ 0.5, the lensing 
efficiency for the model with a constant w - -0.8 is slightly 
smaller than in the ACDM model because the abundance of 
halos is also slightly smaller (see Fig.0}. A similar difference 
appears between the EDE1 and EDE2 models. This is due to 
the fact that in the first the normalisation of the power spec- 
trum is higher than in the second, causing a higher abundance 
of clusters. 

An effect that we also recognise in these plots is that in 
cosmologies with early dark energy, the optical depth per unit 
redshift rises and reaches a significant level already at relatively 
high redshift, while it is still negligible in a ACDM model. As 
discussed before, the models alternative to ACDM that we have 
studied here have a larger fraction of structures at high redshift, 
causing this earlier and larger contribution to the strong-lensing 
efficiency. 

Further detail on this aspect is provided by Fig. [9] In its top 
panels, it shows the cumulative optical depth per unit redshift, 
which we can write as 

C d (z) = P"" X t d (z')dz' , (17) 

normalised to the ACDM case. By its increase towards high 
redshift, it emphasises directly how the lensing efficiency drops 
already at lower redshift in a ACDM Universe with respect to 
the (early) dark-energy cosmologies. The bottom panels show 
the cumulative optical depth per unit redshift normalised to the 
present value in the ACDM model, Cd,ACDM(0) = t^acdm- 
This illustrates the same effect in a different way. For instance, 
we see that the cumulative optical depth per unit redshift in 
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Fig. 8. Similar to Fig.0 but ignoring the effect of cluster merg- 
ers on the lensing efficiency. 

the EDEl model reaches the same value C(1,acdm(0) already 
at z ~ 0.8 that ACDM reaches today. Conversely, the cumula- 
tive optical depth per unit redshift in the ACDM case has al- 
ready dropped by an order of magnitude by z ~ 0.8. Similarly, 
the EDE2 model reaches the total optical depth of the ACDM 
model atz~ 0.7. 

In agreement with our earlier discussion, we note that this 
specific evolution does not depend on whether we take dynam- 
ical processes into account or not. The enhanced lensing ef- 
ficiency in the high redshift tail may have stimulating conse- 
quences, as we shall discuss later. 

The large spikes shown in Fig. are obviously due to the 
variation of the lensing efficiency of galaxy clusters during 
mergers. Very small spikes appear also in Fig. |8j where dy- 
namical processes are not taken into account. There, they orig- 
inate from numerical effects, in particular to the fact that our 
time resolution is very high and the number of haloes is lim- 
ited. Indeed, the spikes become larger well above redshift unity, 
where the number of contributing haloes is reduced (remember 
that each halo is characterised by a different source redshift, 
drawn from a distribution which peaks around z ~ 1.2). 

8. Summary and discussion 

We have analysed the incidence of pronounced (long and thin) 
arcs in galaxy clusters in four dark-energy models. In particu- 
lar, we considered two early dark-energy cosmologies in which 
the density parameter in dark energy at high redshift remains 
small and positive. We compared them to a model with constant 
equation-of-state parameter w de = -0.8 and a ACDM model 
for which Wd e = — 1. 

For each cosmological model, we used Monte-Carlo tech- 
niques to build up merger trees for a set of N — 500 cluster- 
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Fig. 9. Top panels. Logarithm of the cumulative optical depth per unit redshift, normalised to its value in a ACDM universe. 
Curves representing the four cosmologies studied in this work are shown, as labelled in the plot. Bottom panels. Logarithm of 
the cumulative optical depth per unit redshift, normalised to its value at present in a ACDM universe. Left and right panels show 
results including and ignoring cluster mergers, respectively. 



sized dark matter halos. By modelling each halo by an NFW 
density profile with elliptically distorted lensing potential and 
suitably accounting for cluster interactions during mergers, we 
calculated the optical depth per unit redshift both accounting 
for and ignoring cluster mergers. To this end, we also consid- 
ered a realistic distribution for the source redshift. 

We find that, in agreement with the results of iFedeli et alJ 
(2006), cluster mergers enhance the occurrence of arcs by a fac- 
tor between 2 and 3. This occurs in all cosmological models we 
analysed, and the relative increase is approximately the same, 
because the cluster merger rates in the redshift ranges relevant 
for strong lensing (below z ~ 1) are almost identical (see the 
discussion in Sect. 6). 

However, a potentially more important result is that the op- 
tical depth per unit redshift is larger by a factor of ~ 3 in 
early dark-energy models compared to the models with cos- 
mological constant or with a constant equation of state pa- 
rameter Wde = -0.8, while the differences between the lat- 
ter two are close to negligible. There is also a significant dif- 
ference between the two early dark-energy models due to the 
fact that the model EDEl has a higher normalisation parameter 
erg than EDE2 in order to agree with the CMB observations, 
(cf. TableQ. Thus, halos form earlier in model EDEl. This is 
also demonstrated by Figs. |4] and [5] Moreover, the lensing effi- 
ciency drops already at a lower redshift in a ACDM Universe 
than in the different dark-energy models. The optical depth per 



unit redshift has a significant high-redshift tail in early dark- 
energy cosmologies while it is negligible otherwise. 

A main consequence of these results is that they indicate an 
appreciable difference in the incidence of long and thin gravi- 
tational arcs between the ACDM model and models with early 
dark energy. Therefore, arc statistics may provide an interest- 
ing way to investigate into the reliability of these models, al- 
though the precise contribution of Q<ie,sf will probably be bet- 
ter constrained using cluster counts in the X-ray or Sunyaev- 
Zel'dovich regimes, which suffer from lower systematics. 

The presence of early dark-energy, combined with the tran- 
sient boosts due to cluster mergers could help resolve the dis- 
crepancy between the predi cted and observed a bundances of 
gravitational arcs. Since Bartelm ann et alJ Q998I) first pointed 
out this problem for a AC DM universe, much discussion 
developed around this fact ( Meneahetti et 2003b a; 

Wambsg anss et alJl2004t iLi et alJ 120051) . At present, it seems 
that neither the internal structure of the lensing halos nor 
the redshift distribution of the sources can reconcile theory 
and observations. It has been shown here that the effects of 
early dark energy on structure growth interestingly point into 
the right direction. Sim ilar conclusions were drawn also by 
iMeneghetti et alJ J2005t) . where the lensing efficiency of nu- 
merically simulated dark matter haloes in different dark energy 
cosmologies were analysed. Here the haloes are modelled in an 
analytical way, allowing a much higher mass and time resolu- 
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tion. Moreover, the dark energy models studied there were de- 
rived from SUGRA and Ratra-Peebles JPeebles & Ratral20 02) 
potentials, without an early compon ent. In many aspects, ou r 
work is thus complementary to that of Meneghetti et all (12005 ) . 

Finally, the fact that the lensing efficiency in early dark- 
energy models is much higher at high redshift than in the 
ACDM case can be related to the recent unexpected dis- 
covery of high incidence of giant arcs in high reds hift clus- 
ters (Gladders et al.ll2003l IZaritskv & Gonz alez 2003). Future 
searches for strong lensing in distant galaxy clusters may be 
promising to distinguish between cosmological models other 
than the standard ACDM, or at least to gain a deeper under- 
standing of the role of early dark energy. 
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